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Th's  report  represents  the  combined  results  of  the  authors  individual  investigations.  Dr.  Papa 
formulated  the  theory  and  analyzed  the  results  (Section  1.  2.  3,  and  6)  while  Mrs.  Woodworth  developed 
the  mathematical  techniques  (described  in  Section  4)  that  enabled  the  numerical  results  to  be 
obtained  (Section  5). 
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The  Numerical  Evaluation  of  a 
Physical  Optics  Normalized  Cross 
Section  for  a  Rough  Surface 


1.  INTRODUCTION 

This  report  is  devoted  to  two  main  topics:  The  detailed  explanation  of  the  numerical  evaluation  of 
the  physical  optics  normalized  cross  section  of  a  rough  surface  (o°)  and  the  comparison  of  these 
numerical  results  with  two  other  approximate  evaluations  of  the  physical  optics  integral 
representation  for  o°.  In  a  previous  paper1,  the  general  integral  expression  for  the  normalized  cross 
section  a°  was  given  for  a  surface  which  is  randomly  rough  in  one  dimension  only  and  some 
preliminary  comparisons  of  the  three  representations  for  o°  were  made  as  a  function  of  surface 
roughness  and  mean  surface  slope.  In  this  report,  we  analyze  this  general  integral  expression  for  a° 
and  show  how  it  may  be  evaluated  numerically.  Additional  comparisons  are  included. 

The  basic  condition  for  a  physical  optics  (PO)  model  to  be  valid  is  that  R(  »  X  where  Rc  is  the 
average  radius  of  curvature  of  the  surface  and  X  is  the  em  wavelength.  In  previous  work2,  we  showed 
that  the  condition  T  »  X  (T  is  the  surface  correlation  length)  is  a  sufficient  condition  for  Rc  »  X  if  the 
surface  slopes  are  small,  but  it  becomes  necessary  as  well  as  sufficient  if  larger  slopes  are  not  excluded. 
That  demonstration  assumed  Gaussian  heights  and  correlation.  In  this  paper,  we  will  be  concerned 


(Received  for  publication  16  February  1988) 

1.  Papa.  R.J.  and  Lennon,  J.F.  (1988)  Regions  of  Validity  for  Some  Rough  Surface  Scattering 
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mainly  with  small  slope  conditions.  Under  the  assumption  Rc  »  X.  Barrtck  and  Peake3  have  given  a 
general  expression  for  a  PO  o°  with  no  shadowing.  This  expression  generally4 5 6*7  has  been  evaluated  by 
taking  the  high  frequency  limit,  k  =  2x/X  -»  °°  (GO)  which  allows  decorrelation  of  the  ensemble 
averaging  over  the  height  distribution  from  that  of  the  surface  slope  distribution  followed  by  an 
asymptotic  algebraic  result.  It  has  been  shown1  that  their  product  form  representation: 

<r°  =  I  l2J  (1) 

is  valid  without  the  restriction  X  -*  O.  so  long  as  the  surface  slopes  are  small  (o/T  <  1)  and  J  is  kept  as 
an  integral  form.  Here.  (5^  is  the  scattering  matrix  element  (8),  p  refers  to  the  polarization  of  the 
incident  wave  and  q  refers  to  the  polarization  of  the  scattered  wave  (horizontal  or  vertical).  J  is 
proportional  to  the  probability  density  function  of  the  surface  slopes  and  a  is  the  standard  deviation 
in  surface  height.  Gaussian  heights  and  correlation  are  again  assumed.  This  result  and  the  high 
frequency  form  are  then  compared  with  numerical  evaluation  of  the  general  expression  for  o°  for 
forward  scattering  and  horizontal  polarization. 

The  conventional  expression  for  the  normalized  cross  section  o°  may  be  defined  by  applying 
physical  optics  principles  to  rough  surface  scattering,  where  the  Klrchhoff  integral  is  used  to  represent 
the  scattered  em  wave  and  the  boundary  conditions  on  the  surface  are  satisfied  by  employing  the 
Fresnel  plane  wave  reflection  coefficients1'8.  Multiple  scattering  is  not  Included.  Following  these 
procedures  the  Barrick  and  Peake3  generalized  expression  for  the  field  scattered  from  a  rough  surface 
is  (horizontal  polarization): 

,  L/2  L/2 

glnc  _  -ik|elkRo/4xR0|Eh  J  j  F(4*.  5y)ie,k|k'-ksl  r]dxdy  (2) 

-L/2  -L/2 


where 

Eh  =  incident,  horizotally  polarized  electric  field. 
R0  =  distance  from  origin  to  observation  point. 


3.  Barrick.  D.E.  and  Peake.  W.H.  (1967)  Scattering  from  Surfaces  with  Different  Roughness  Scales: 

Analysis  and  Interpretation,  Battelle  Report,  AD662751. 

4.  Semyonov.  B.I.  (1966)  Approximate  computation  of  scattering  of  electromagnetic  waves  by 

rough  surface  contours,  Radiotekhnika  lElectrontka  11,  1179-1187. 

5.  Barrick.  D.E.  ( 1968)  Rough  surface  scattering  based  on  the  specular  point  theory,  IEEE  Trans. 

Antennas  Propag.  AP-16(4),  449-454. 

6.  Kodis,  R.  ( 1966)  A  note  on  the  theory  of  scattering  from  an  irregular  surface.  IEEE  Trans. 

Antennas  Propag.  AP-14(1),  77-82. 

7.  Sancer,  M.I.  (1969)  Shadow  corrected  electromagnetic  scattering  from  a  randomly  rough  surface, 

IEEE  Trans.  Antennas  Propag.  AP-17(5),  575-585. 

8.  Ruck.  G.T..  Barrick,  D.E..  Stuart.  W.D.,  and  Krichbaum.  C.K.  (1970)  Radar  Cross  Section 

Handbook,  VoL  2,  New  York.  Plenum  Press. 
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r  =  xx  +  yy  +  4(x,y)s  =  distance  from  origin  to  point  on  rough  surface. 


kt,  k8  =  unit  constant  vectors  pointing  in  direction  of  Incidence  and  scattering. 

4X.  4y  =  local  surface  slopes  in  x  and  y  directions  at  surface  point 
4(x,y),  that  is,  34 /3x  and  34 /3y . 

The  factor  F(  4X.  )  is  a  function  of  the  local  normal  to  the  surface  and  of  local  Fresnel  reflection 

coefficients  at  each  surface  point.  The  Barrick  and  Peake  expression  for  F  is  incorrect  as  printed.  The 
correct  form  is  given  by  Sancer7. 

To  form  the  normalized  cross  section  of  the  rough  surface.  o°.  it  is  necessary  to  calculate 

<y°  a  <ES*  •  Es>  -  <ES>2  (3) 

where  <  •  >  denotes  an  ensemble  average  over  the  random  variables  4-  4i,  4x-  4y  4xi .  and  qyl .  Here,  4  is 
the  random  hight  at  point  (x,  y)  and  4i  is  the  random  height  at  point  (x, .  y j).  The  general  PO  expression 

(Eq.  (3)1  for  o°  involves  a  ten-fold  integration  over  the  variables  x.  y,  x,.  y,.  4.  4i.  4X' 4y,  4xi  •  ar“d  4yi-  By 

making  use  of  the  stationarity  of  the  random  process,  this  expression  can  be  reduced  to  an  eight-fold 
integral.  By  what  is  traditionally  referred  to  as  a  geometrical  optics  assumption,  this  result  is  further 
reduced  to  a  single  integral  form  and  finally  evaluated  as  an  algebraic  expression..  The  single  integral 
form  can  also  be  obtained  by  making  use  of  a  small  slope  condition1.  For  that  case,  though,  the  final 
evaluation  does  not  apply. 


2.  THEORY  FOR  SCATTERING  FROM  A  ONE-DIMENSIONAL  ROUGH  SURFACE 

To  simplify  the  general  expression  for  o°  and  to  make  numerical  evaluations  of  a°  more  efficient, 
it  will  be  assumed  that  the  random  distribution  in  heights.  4.  have  only  a  one-dimensional  variation 
4  =  4M-  It  will  also  be  assumed  that  the  scattering  takes  place  in  the  direct  forward  direction,  so  that 
there  is  no  azimuthal  vaiiaaoi.  ivs  -  0°)  and  that  there  Is  shadowing  The  most  peneral  expression 
for  o°  for  a  one-dimensional  rough  surface  involves  a  six-fold  integral,  where  the  variables  of 
Integration  are  x,  Xj.  4p  42-  hi-  and  Here.  X]  and  Xj  ,  are  two  points  on  the  rough  surface.  4i  and  42  are 
the  random  surface  heights  at  the  two  points  and  p,  and  Pj  are  the  random  surface  slopes  at  the  two 
points.  By  using  the  fact  that  the  surface  heights  are  to  be  regarded  as  a  stationary  random  process,  so 
that  the  correlation  function  of  the  surface  heights.  c(t1.  is  a  function  only  of  the  separation  between 
points  x  =  xj  -  X2  .  the  six-fold  integral  for  o°  may  be  reduced  to  a  five-fold  integral.  Then,  using  the  iact 
that  the  general  expression  for  cr°  is  a  function  only  of  the  height  differences  4  =  4i  -  42'  the  five-fold 
Integral  can  be  reduced  to  a  four-fold  integral.  Finally,  by  assuming  the  trtvanate  distribution 
function  for  the  heights  4  and  slopes  P]  and  gj  is  a  Gaussian  function,  the  integration  over  4  can  be 
performed  analytically.  Then,  the  remaining  expression  for  o°  reduces  to  a  triple  integral,  which  can 
be  performed  numerically  using  Gaussian  quadrature  techniques. 
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Following  the  development  given  by  Hagfors9,  the  expression  for  a°  may  be  written  as  (1): 

oo 

c°  =  2itk  Jdx  cosvKx  w(x)  (4) 


where 

vx  =  klslnO,  -  sin6s) 

r  =  x  1  -  x2  (separation  between  two  points  on  the  rough  surface) 

and 

w(t)  =  G(x)  -  H.  (5) 


G(x)  =  2 


oo  oo  oo 

Jdp2F(p2)  Jdp,F(pi)  Jd^  cosvz  %  •  P3(p,.  p2.  s-  *) 

—  oo  —  oo  0 


P3  =  trivariate  distribution  function  in  height  differences  and  surface  slopes.  In  addition,  s  i  .  M  i  • 
and  ^2*  h2  are  the  heights  and  slopes  of  point  1  and  point  2  .  respectively  ,  with  %  -  t, ,  -  s2.  M  i  =  s*  i  and 
1*2  =  ^x2- 

F(p)  is  a  complicated  function  of  the  slopes  p.  and  is  given  in  the  appendix,  as  is  the  form  for  the 
quantity  H.  It  can  be  shown  that  w(x)  -*  0  as  x  —  <*>. 

There  are  no  restrictions  on  the  expression  for  0°  given  by  |Eq.  (4)1.  other  than  the  validity  of 
physical  optics.  T  »  X.  There  are  singularities  in  the  integrand  of  )Eq.  (4))  but  they  can  be  shown  lo  be 
integrable  and  the  expression  has  been  evaluated  numerically,  using  quadrature  formulas  to  provide  a 
standard  for  the  two  limiting  case  solutions.  For  the  two  sets  of  assumptions  considered  here,  the 
expressions  tn  (Eq.  (4)1  and  (Eq.  (5)1  reduced  to  the  simpler  form  oflEq.  (6)]  involving  a  single 
integration  over  x: 

oo 

o°  =  2itk  F2(uap)  Jdx  cosvxx  •  1*2  -  X?l  (6) 

o 

where 

X2  =  expi-c^zll  -  p)i 


p  =  exp(-x2/T2) 

X  i  =  exp(-v2(J  2/2) 


9.  Hagfors.  T.  (1964)  Backscatlering  from  an  undulating  surface  with  applications  to  radar  returns 
from  the  moon. Geophysical  Res.  69(18),  3779-3784. 
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and 


v =  -k(cos0,  +  r-;  *s). 

The  details  of  the  arguments  In  each  case  are  presented  in  the  appendix.  Finally,  for  the  case  where 
X  -»  0  and  L  ->  1.  the  integral  of  (Eq.  (6)|  can  be  explicitly  evaluated  for  a0  using  the  stationary  phase 
methed'’: 

cr°  =  7tkF2(psp)(T  ViT/Slexpl  -v^T2/4I2)  (7) 

Here.  L  =  -  ovz 

One  of  our  objectives  here  is  to  examine  the  conditions  under  which  each  of  these  two  forms  agrees 
with  a  direct  numerical  evaluation  of  (Eq.  (4]|  and  (Eq.  (5)1  and  to  define  the  regions  where  one  or  the 
other  limiting  case  solutions  would  be  preferred.  These  are  discussed  in  Section  5. 


3.  SIMPLIFICATION  OF  EXPRESSION  FOR  NORMALIZED  CROSS  SECTION 


|Eqs.  (4)  and  (5)1  constitute  exact  physical  optics  expressions  for  the  normalized  cross  section  a  of 
a  one  dimensional  rough  surface  These  equations  show  that  this  general  expression  for  involves  a 
four-fold  integration  over  the  variables  r.  4-  p,  and  |i2-  ^  l^c  Irivariale  distribution  function  P3fij.P1. 
p2)  in  4.  p  |  and  p2  is  assumed  to  be  a  Gaussian,  then  the  integration  over  4  may  be  accomplished 
analytically,  as  shown  in  the  appendix.  (Eq.  A5).  that  is. 


=  |d4cosv,4  •  P<  =  —  , - 

0J  (2nVM„) 


cosv  Be  V,2/4A 


(8) 


where 

A  =  M,,/(2IRI) 


B  =  M)2/M|  i  )(p  |  +  p2) 

c  =  (1/2IRI)  UM22  -  M22/Mu)(pi  +  p 2}  +  2(M 23  -  Mi2/Mn)p,p2l 

is  the  co-factor  of  the  covariance  matrix  R  tj.  One  may  rewrite  the  quantity  c  as 

c  =  (a2/T2)(l/Mu)(E(pf  +  p2)  +  2FpiP2| 


where 

E  =  (l/(a2/T2)  {(MuM22  -  M?21/(2IRI)) 


and 
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F  =  (1/(C2/T2)  ((MnM23  -  M?2  )  /(2  I R I )) 


After  considerable  algebraic  manipulation,  one  can  show  that  E  =  1  and 
F  =  -exp(-x2/T2)(  1  -  2t2/T4) 

1  -  F2  =  Mu/(4ct4/T4) 


By  completing  the  square,  one  may  rewrite  c  as 


c  =  (a2/T2)(  1/Mu) 


-  2, 


(hi  +  42  +  2Fhih2  +  ^hl) 


(9) 


Now.  define  =  h2  +  Fhj.  so  that  c  becomes 

<j2  p22  ^ 

C  "  L(T2)(Mu  1  +  4(cr2/T:i)  . 

The  question  arises,  what  happens  to  the  expression  for  I  given  by  (Eq.  (8)1  as  t  -»  0  (x1 

24(T^  _  ftryS 

can  be  shown  that  Mn -♦  (-^nj—)(x2/T2)  and  IRI  t_"0  )  (t8/T8).  Also,  as  t  -*  0. 

0.  F—  -1  and  pt  -•  so  that  p’2  0.  Then  we  have 

c  ^T0  [c,  +  hi/4(o2/4(a2/T2))]. 

where  c ,  is  a  constant. 

So  that  if 

hi  =  0,  c  tZTQ  c, 


(10) 
x2).  It 


and 

I  rr0  (l/(2x  Ir  I )e~c l 

Hence.  I  and  therefore  w(x)  has  a  singularity  at  x  =  0.  However,  it  can  be  shown  that  this  singularity  is 
tntegrable.  This  may  be  seen  from  (Eq.  (5)1.  Here,  the  behavior  of  W(x)  near  x  =  0  is  governed  by  G(x). 
which  in  turn,  is  governed  by  the  behavior  of  P3  near  x  =  0.  But,  from  [Eq.  (A3)!  in  the  appendix. 

P3  =  (2tx) -3/2  I R  I -l/2e-© 
where 

9  =  uT  R-‘u. 

Now. 

IRI  (8<J6/T4)(x8/T 8) 


and 
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9  t-r~o  l/(c2x8)[c3x2£2  +  c4x34(h1  +  p2). 


where  c2  c3,  and  c4  are  constants. 

Therefore,  since  %  —  0  and  hi  9  T“Z*o  +  •» 
so  that 

P3  t  ->’o  0 

and  W(t)  is  llnlnte  as  x  -» 0.  Thus.  W(x)  is  integrable  about  the  point  x  =  0.  How  this  may  be 
accomplished  numerically  is  explained  in  Section  4. 
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a  =  (1  -  x2/T2  +  V2x4/T4  -e-*2^2)  tTT0  -V6t6/T6 

P  =  (1  +  x2/T2  +  V2x4/T4  -e«2/T2)  X-ZT0  V6x6/T6 

Note  the  first  three  elements  of  a  and  3  are  the  first  three  terms  of  the  Taylor  series  expansion  of  the 
x2  x2 

exponential.  Let  x  =  -  for  a  and  x  =  +  for  3-  The  Taylor  series  expression  for  ex  Is 

a 

Zxh 

hT 

h*0 


50  = 


h=3 


if  Ixl  is  made  small  ex  =  1.  For  example  if  x  =  0.01  then 
10~4  10“6 

e-°-01  =  1  -  0.01+  — ig - -  +  .  .  . 

Many  or  all  of  the  significant  figures  could  be  lost  by  finding  ex  then  subtracting  the  first  few  terms. 
So  the  program  was  written  to  find  (1  -  ex  )  and  ( 1  +  x  +  x2/2  -  ex)  directly  by  the  sum  series  for 

oo 

V1  Xh 

small  Ixl.  2_j  "JJT  was  used  for  M  u  where 
h=l 


M, 


=  4?pr  £(  1  -  e'2T2/T2)  +  4e'2T2/T2(  -  =fT )]  r^o  24  fr(z2/T2) 


(12) 


4. 1  Limits  of  Integration 

The  function  drops  off  exponentially  in  and  p2  so  we  use  this  property  to  determine  the 
practical  limits  of  integration.  First  we  integrate  [Eq.  (A  1)1  to  get  H.  The  exponential  part  is 

exp  -  [p/i2°/x)]2 

The  function  drops  off  by  e~16  =  10~7  when  =  +  8a/x 
As 

t  -*  «o.  G(x)  -*  H  so  w(x)  =  G(x)  -  H  -»  0 

In  some  cases  H  is  extremely  small  so  it  can  be  neglected  and  the  behavior  of  only  G  need  be  considered. 
For  single  integral  physical  optics  (PO)  we  use  an  upper  limit  of  zmax  =  4T.  At  that  point  the  difference 
has  dropped  off  by  e-16  ■»  lO-7. 
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f  p  *  '  P"  ' » 


The  full  triple  Integral  for  physical  optics  (TI)  is  a  very  complex  function  in  t.  Therefore  the 
limits  cannot  be  determined  analytically.  However,  we  do  expect  a  limit  comparable  to  the  PO 
integration  limit.  We  find  the  value  of  w(x)  for 

T 

x  =  n  y  for  n  =  1.2 . 20 

The  first  point  where  two  consecutive  values  of  w(x)  drop  below  10~io  is  taken  as  .  the  upper  limit 
of  integration.  Two  consecutive  w(xn)  are  checked  in  case  one  is  a  null  caused  by  the  function's 
oscillations.  Our  experience  showed  this  method  produces  reasonable  values  for  xnvu(  . 

Next  we  find  the  limits  in  p,  and  p2  .  The  exponential  part  of  the  expression  is  e~r  where 

c  =  +  (  p,/(2c/x))2 

Good  limits  are 

Pimax  =  ±  a/r 

p-2  =  ±  4ifM^/(o/T) 

Either  will  cause  c  to  Increase  by  16  which  causes  the  function  to  drop  olTby  e  16  =  10  7.  Note  that 
p2max  changes  with  x  sinced  M, ,  is  a  function  of  x.  When  x  is  small  M , ,  and  p2n)_.lx  are  small  so  the 
function  is  integrated  over  a  narrow  range.  When  x  is  large  the  range  of  integration  in  p'2  is  broader. 
The  flexible  integration  range  cuts  down  on  wasted  computational  effort  at  small  x. 

4.2  Integration  Procedure 

The  next  step  is  to  pick  an  integration  scheme.  Common  integration  schemes  are  Simpson's  rule 
and  Romberg  integration  which  is  a  higher  order  variation  on  Simpson's  rule.  Unfortunately,  the 
function  oscillates  and  the  results  diverge  even  with  more  than  16.000  points  and  days  of  CPU.  This 
occurs  because  Simpson's  rule  acts  like  a  digital  filter  amplifying  the  large  high  order  frequencies.  At 
finer  step  sizes,  the  Integrated  sum  will  grow  in  size  for  a  high-frequency  function  (Hamming10  .  p.  39). 

Another  common  integration  scheme  is  Gaussian  quadrature.  Plots  of  the  function  in  p,  and 
p2  show  that  it  has  about  10  to  20  cycles  of  oscillation  in  the  Interval  of  integration.  So  each 
integrand  was  divided  up  into  20  sections  and  each  section  was  solved  by  Gauss-Legendre  quadrature 
and  the  sections  were  totaled  up  to  give  the  integral  for  that  integrand. 

The  integration  is  done  in  a  very  complex  iterative  method.  Before  integrating,  the  abscissas  and 
weights  are  determined  for  Gauss-Legendre  integration  with  2.  4.  8.  16.  32.  and  64  points.  Each  section 
was  first  found  by  2  point  quadrature,  then  with  4  point,  then  8  point,  etc.  After  each  step  the  sections 
and  their  sum  are  compared  to  see  if  the  results  have  converged.  The  sections  which  converge  or  fall 
below  the  error  threshold  are  marked  "done".  When  the  total  of  the  sections  converges  we  have  the 
solution  for  that  integral.  This  is  done  for  all  three  levels  of  integration. 


10.  Hamming,  R.W.  (1977)  Digital  Fillers,  Prentice  Hall.  New  Jersey. 
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In  the  outer  two  integrations  (over  r  and  Hi),  an  error  threshold  is  determined  which  is  passed  lo 
the  next  inner  integration  (over  p.l  and  p'2  ).  If  a  section  or  the  sum  falls  below  this  threshold,  the 
respective  portion  of  the  inner  integration  is  marked  "done"  even  if  it  has  not  converged,  because  the 
value  is  too  small  to  significantly  contribute  to  the  final  result.  The  complicated  iterative  method 
saves  computer  time,  but  greatly  increases  the  program  complexity. 

Table  1  shows  the  subroutine  calls  needed  to  perform  the  actual  triple  integration.  The  program  is 
ten  subroutines  deep.  MAIN  is  the  main  program.  It  calls  subroutine  QUADT.  Subroutines  QUADT, 
QUADM1.  and  QUADM2  set  up  the  parameters  for  the  quadrature  over  x.  p[  and  p2  respectively.  The 
quad  subroutines  call  subroutines  1G,  1GM1.  and  1GM2  respectively.  These  subroutines  do  the  actual 
Integrations  of  each  section  of  the  integrand  over  x.  plt  and  p2.  The  1G  subroutines  respectively  call 
wcos,  R.  and  P.  which  find  the  functions  in  x.  and  p2 .  wcos  and  R  call  the  next  inner-most 
quadrature  subroutines  (QUADM1  and  QUADM2).  Subroutine  WCOS  also  calls  Si  and  S3  which  find 
( 1  -  e*)  and  ( 1  +  x  +  x2/2  -  e*)  respectively.  Subroutines  R  and  P  call  F  which  finds  F(p). 

-im  Table  1.  Calculation  of  Normalized  Cross  Section  o° 

Geometric  Optics 


5.  NUMERICAL  RESULTS 

In  this  section,  graphs  of  o°  will  be  presented  as  a  function  of  the  Rayleigh  roughness  parameter  I 
for  different  slope  conditions,  o/T.  For  a  one  dimensional  rough  surface,  three  models  are  used  to 
represent  the  normalized  cross  section,  o':  (1)  The  triple  Integral  representation  (TI)  given  by  (Eqs.  (4). 
(5),  A1  and  A5.)  (2)  the  single  integral  representation  (PO)  given  by  (Eq.  (6)).  valid  for  small  slopes  and 
(3)  the  asymptotic,  high  frequency  representation  (GO),  given  by  (Eq.  (7)).  In  this  report,  the  parameters 
used  to  generate  the  figures  were  taken  to  have  the  following  values:  \  ranges  from  0.01  m  to  0.35  m. 


MAIN 

I 

F 


Physical  Optics 

Triple  Integral 

MAIN 

1 

MAIN 

1 

1 

QUAD 

1 

1 

QUADT 

1 

IG 

1 

1 

IG 

1 

1 

FNPO 

1 

S3— WCOS— Si 

1 

1 

QUADM1 

1 

1 

IGM1 

1 

1 

R-F 

1 

1 

QUADM2 

1 

1 

IGM2 

4>g  =  0°  and  the  complex  dielectric  constant  of  the  surface  e  =  4.0.  The  Rayleigh  parameter  was  varied 
by  varying  the  em  wave-length  X.  The  angles  of  Incidence  and  scattering  are  always  fixed  at  89.75°. 

In  Figure  1.  the  mean  slope  a/T  =  0. 1,  so  that  the  PO  single  integral  representation  is  valid.  It  can 
be  seen  that  the  PO  representation  of  <f  is  within  0. 1  percent  of  the  exact  TI  representation  for  all 
Rayleigh  parameters  X.  The  high  frequency  GO  representation  is  within  0.5  percent  of  the  TI 
representation  only  for  X  >  5.  For  I  <  2,  the  GO  representation  ts  completely  inaccurate. 


SLOPE  0.1 


I 


Figure  1 .  Normalized  Cross  Section  o°  vs.  Rayleigh  Parameter  X.  a/T  =  0.1 
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In  Figure  2,  Che  mean  slope  cr/T  has  been  increased  to  0.2.  Now.  there  is  a  greater  discrepancy 
between  the  PO  representation  and  theTl  representation,  but  the  PO  representation  ts  still  within  0.5 
percent  of  the  Tl  representation  for  all  values  of  I.  The  GO  representation  is  inaccurate  for  £  <  4. 
compared  to  the  Tl  representation. 


SLOPE  0.2 


I 


Figure  2.  Normalized  Cross  Section  vs.  Rayleigh  Parameter  £.  c/T=  0.2 


12 


In  Figure  3.  the  mean  slope  has  been  increased  to  0.5.  Here,  the  discrepancy  between  the  PO 
representation  and  the  TI  has  increased  to  about  5  percent  for  I  >  5  and  is  as  large  as  20  percent  in  the 
vicinity  of  I «  2.  One  can  infer  that  the  single  integral  PO  representation  is  accurate  only  for  slopes 
o/T  <  0.5.  The  GO  representation  is  within  5  percent  of  the  accurate  Tl  representation  only  for  £  >  5. 


SLOPE  0.5 


I 


Figure  3.  Normalized  Cross  Section  o°  vs.  Rayleigh  Parameter  I.  a/T  =  0.5 
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In  Figure  4,  the  mean  slope  has  Increased  to  1.0.  The  discrepancy  between  the  PO  representation 
and  the  TI  is  20  percent  for  I  >  5  and  is  as  large  as  45  percent  near  1=3.  Th  s,  the  single  integral  PO 
representation  is  not  valid  or  accurate  for  c/T  >  0.5.  Also,  It  should  be  noted  that  the  high  frequency  GO 
representation  is  not  accurate  for  slopes  a/T  >  0.5,  regardless  of  the  value  of  the  Rayleigh  parameter  I. 


SLOPE  1 


Figure  4.  Normalized  Cross  Section  o°  vs.  Rayleigh  Parameter  X,  o/T  =1.0 


In  Figure  5,  the  main  slope  has  been  increased  to  2.0.  Here,  the  discrepancy  between  the  single 
Integral  PO  model  and  the  TI  is  100  percent  for  a  range  of  t  values,  so  that  neither  the  PO  model  nor  the 
GO  model  are  accurate  for  most  I  values. 


SLOPE  2 


I 


Figure  5.  Normalized  Cross  Section  a°  vs.  Rayleigh  Parameter  I,  o/T  =  2.0 


T 


In  Figure  7.  plots  are  shown  of  A  (In  percent)  vs  o/T  for  different  Rayleigh  parameter  £  regimes, 
that  Is,  for  different  frequency  ranges.  Here,  the  percent  error 

o°(TI)  -  o°(PO) 

A  “  o°(PO) 

For  I  <  1  we  took  £-  1/2,  for  1  <  £  <  4  we  took  I »  3,  and  for  £  >  4  we  took  1-6.  It  may  be  noted  that 
the  greatest  error  lies  In  the  Rayleigh  parameter  regime  1  <  £  <  4.  Also,  the  error  increases  as  the  slope 
o/T  Increases.  If  one  calculates  according  to  the  single  Integral  (PO)  representation,  one  can  recover 
the  exact  TI  o°  value  by  multiplying  the  cr°(PO)  value  by  the  appropriate  A  value  and  adding  it  to  the 
(PO)  value,  that  Is 

o°(TI)  =  o°(PO)  +  A|c°(PO)]. 


ERROR  OF  PO  APPROXIMATION 


Figure  7.  Relative  Error  vs.  Slope  for  Different  Rayleigh  Parameters 
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&  DISCUSSION 


In  this  report,  it  was  shown  that  an  exact  representation  of  a  physical  optics  model  for  the 
normalized  cross  section  c°  could  be  derived  for  a  one-dimensional  rough  surface.  This  exact 
representation  of  <r°  was  reduced  to  a  triple  integral  over  the  separation  t  and  the  surface  slopes  at  two 
points.  Pi  and  g2.  It  was  demonstrated  that  the  integral  over  x  had  a  singularity,  but  that  it  was 
integrable.  The  integrals  could  be  performed  numerically.  The  Romberg  method  of  numerical 
integration  did  not  converge;  a  Gaussian  quadrature  technique  had  to  be  employed.  Three  models  for 
a°,  the  triple  integral  (TI).  single  integral  physical  optics  (PO).  and  the  geometrical  optics  (GO)  were 
studied  as  a  function  of  Rayleigh  parameter  I  for  different  slope  (a/T)  regimes.  The  Rayleigh 
parameter  was  varied  by  varying  the  frequency;  the  angles  of  incidence  8,  and  scattering  6S  are  fixed  at 
89.75°.  A  graph  of  relative  error  A  vs  <r/T  for  different  Rayleigh  parameter  I  regimes  was  presented.  It 
can  be  used  to  recover  the  exact  o°  (TI)  value  from  the  single  integral  PO  representation  of  o°,  for  this 
set  of  conditions. 


7.  CONCLUSIONS 

It  was  found  that  the  PO  model  is  accurate  only  for  a/T  <  0.5.  regardless  of  the  Rayleigh  parameter 
It  was  also  found  that  taking  the  high  frequency  limit  X-*  0  is  not  a  sufficient  condition  for  the 
validity  of  the  GO  model.  From  the  numerical  results,  it  was  shown  that  the  surface  slopes  must  be 
small;  a/T  <  0.5.  From  analysis  of  the  integral  representation  of  a°.  it  can  be  shown  that  the 
asymptotic  expression  (GO)  for  a°  may  be  derived  if  I  >  4. 
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Appendix 

The  Relationships  in  the  General  Integral  Expression  for  Cross  Section 


The  purpose  of  this  appendix  is  to  give  the  details  of  the  respective  arguments  by  which  the  results 
of  [Eq.  (4)|  and  [Eq.  (5)1  can  be  simplified  to  those  of  (Eq.  (6)1  for  both  sets  of  assumptions.  The  first  step 
Is  to  examine  the  elements  of  the  integrands  in  more  detail. 

First,  the  quantity  H  used  in  w(r)  In  |Eq.  (4)J  Is  given  by  the  expression 

r  °o  ■■2 

t2, 


H  = 


rT2e-Vj,°2/(4xa2)l  fexp  -  (p/(2a/T)  )2F(p)dp 

*  /  _  OO 


(Al) 


The  function  F(p)  used  tn  [Eq.  (5)1  can  be  reduced  to  the  form 

F(p)  =  (1  -  R  HhsinO,  +  cos9()  +  [1  +  R  (y , ) )  (psln0s  -  cos0s)  IA2) 


where  it  has  been  assumed  that  there  is  no  y  variation  ( r^-=  0 )  and  the  scattering  takes  place  In  the 

dy 

forward  direction  (0S  =  0°).  Here,  the  Fresnel  reflection  coefficient  is  given  by 
R  x<Yi)  =  (cosy,  -  -Je  -sln2y,)/(cosy,  +  Ve  -sin2y,) 


and 

cosy,  =  (psin0,  +  cos0,)/(V  1  +  2'p2). 

The  trivariate  distribution  function  P3  in  [Eq.  (5)]  is  given  by 
P3  =  (2n ) -3/2  iri-i/2  e-5fTi*_I 

tjT  R-,u  is  a  positive  definite  quadratic  form: 


(A3) 


uT  R  !u  =  1/(2  IBHIMj,  +  2M12  +  p2)  +  M22  (p?  +  P-i)  +  2M23p1p2l 


uT  =  (4.  Up  p2) 

I R I  is  the  determinant  of  the  surface  height  covariance  matrix 


f2[o2  -  p) 

R  =  -dp/dt 

-dp/at 


a2p/at2  it  =  0 

-a2p/ar2 


-ap/at 
-a2p/ar2 
+a2p/ar2  it . 0 


where 


p  =  surface  correlation  function 
p  =  <j2exp(-r2/T2) 

It  should  be  noted  that  the  expression  for  R  given  by  Hagfors9  has  errors  in  the  signs  of  several 
elements. 

M  j j  is  the  co-factor  of  the  covariance  matrix  R,  j .  The  triple  integral  in  (Eq.  (5))  may  be  reduced  to 
a  double  integral  when  t  *  0  by  using  the  known  expression  for  the  cosine  transform  of  a  Gaussian 
function 


Jd^cosvz  3  =  1  /(2n  Vm,  i  )e~c  cos(vzB)  expl-vi?/4A) 


where 


A  =  Mi.  /  ( 2  I  R  I ) 


B  —  (Mj2/Mjj)(p|  +  p2) 

C  =  (1/2IRI)  J(M22  -  M?2/M ,, )(m?  +  P2)  +  2(M23  -  M^2/Mm)PiP2J 

At  this  point  we  turn  to  the  two  cases  that  are  derived  under  different  assumptions. 

First,  consider  the  small  slope  case,  o/T«l .  In  this  case,  the  covariance  matrix  R  becomes 

f  2(o2  -  p|  0  0  \ 

R  =  0  0  0  (A6: 

^  0  0  0  , 

so  (hat  the  heights  and  slopes  are  decorrelated.  Then  the  integrations  over  p,  am  p2  111 IEQ-  (5)1  may  be 
approximated  by  setting  up  pj  =  p2  =  constant  =  tan  y  =  psp  where  tan  y  is  the  slope  of  a  facet  which 
will  produce  a  reflected  wave  in  the  specular  direction: 
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tany  =  |  sln0,  -  sln0s  |  /(cos,  +  iosfls)|A7| 

(Eq.  (4)1  for  or°  now  reduces  to  a  single  integration  Ibnn  as  shown  In  |Eq.  (fill.  It  should  he  noted  that 
since  the  slopes  are  assumed  to  be  small  (a/T  «.  1).  then  p  =  tan  y  so  that  y  ^  0.  which  implies  the 
specular  condition  08  =  0,. 

The  argument  for  reducing  the  integrals  In  |Eq.  (5)]  to  the  result  in  |Eq.  (6)1  is  somewhat  different 
for  the  high  frequency  case.  In  the  GO  limit  the  integrals  over  p ,  and  p2  maV  be  approximated  bv 
removing  F(pj)  and  F(p2)  from  the  integrals  and  setting  them  equal  to  constants;  the  justification  for 
this  is  the  stationary  phase  (or  specular  point)  argument.  This  argument  states  (hat  for  large  k  in  the 
exponential  (or  cosine)  factor,  the  only  surface  regions  which  contribute  to  the  integral  are  those 
smoothly  curving  portions  in  a  position  to  specularly  reflect  into  the  desired  scattering  direction. 
Then,  pt  »  =  tan  Tf  =  Psp  and  lEqs.  (4)  and  (5)1  reduce  to  [Eli]  (6)|  (see  Barrick  and  Peake  (3)). 
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MISSION 

Rome  Air  Development  Center 

RAVC  plam  and  exec utes  research,  development,  test 
and  selected  acquisition  programs  In  support  o< 
Command,  Conteot,  Communications  and  Intelligence 
( C  * T )  activities .  Technical  and  engineering 
support  within  areas  o{  competence  Is  provided  to 
ESP  Program  OUlces  (P0*)  and  otfcea  ESP  element* 
to  perform  elective  acquisition  ot  CSI  systems. 

The  areas  o£  technical  competence  Include 
communications ,  command  and  c ontrol,  battle 
management,  information  processing,  surveillance 
sensors,  Intelligence  data  collection  and  handling, 
solid  state  sciences,  eEecteomaonet<cc4,  and 
propagation,  and  electronic,  maintainability, 
and  compatibility. 


